source("C:\\projects\\quantileMap\\R code\\release\\quantilemap.R")

set.seed(12345)
ndist<-400
xdim<-20
mu<-0.5
sigma<-0.1
lambda<-10
n<-rpois(ndist, lambda)
data<-as.list(rep(NA, ndist))
for(i in 1:length(data))
{
	if(n[i]>0)
		data[[i]]<-rnorm(n[i], mu, sqrt(sigma))
}

data.EB<-qm.EB.FI(data)
data.MLE<-qm.MLE.FI(data)

qm.plot(data, xdim, FI=FI.EB, col.palette="RGB")
qm.plot(data, xdim, FI=FI.MLE, col.palette="HCL")
